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Ground state cooling of massive mechanical objects remains a difficult task restricted by the unresolved me¬ 
chanical sidebands. We propose an optomechanically-induced-transparency cooling scheme to achieve ground 
state cooling of mechanical motion without the resolved sideband condition in a pure optomechanical system 
with two mechanical modes coupled to the same optical cavity mode. We show that ground state cooling is 
achievable for sideband resolution oJm/K as low as 0.003. This provides a new route for quantum manipula¬ 
tion of massive macroscopic devices and high-precision measurements. 
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INTRODUCTION 


Cavity optomechanics provides a perfect platform not only 
for the fundamental study of quantum theory but also for the 
broad applications in quantum information processing and 
high-precision metrology [1-3]. For most applications it is 
highly desirable to cool the mechanical motion to the quan¬ 
tum ground state by suppressing thermal noise. In the past 
few years numerous efforts have made strides towards this 
goal through backaction cooling [4-13]. However, the cool¬ 
ing limit is subjected to quantum backaction, and ground state 
cooling is possible only in the resolved sideband (good-cavity) 
limit [14, 15], which requires the resonance frequency of the 
mechanical motion (Wm) to be larger than the cavity decay rate 
K. This sets a major obstacle for the ground state preparation 
and quantum manipulation of macroscopic and mesoscopic 
mechanical resonators with typically low mechanical reso¬ 
nance frequency. Therefore, it is essential to overcome this 
limitation, so that ground state cooling can be achieved irre¬ 
spective of mechanical resonance frequency and cavity damp¬ 
ing. 

Some recent proposals [16] focus on circumventing the 
resolved sideband restriction by using dissipative coupling 
mechanism [17, 18], parameter modulations [19-23] and 
hybrid systems [24-29]. Here we propose a unresolved- 
sideband ground-state cooling scheme in a generic optome¬ 
chanical system which does not require modified mechanisms 
of coupling or specific modulation of the system parameters, 
or additional components. We take advantage of the destruc¬ 
tive quantum interference in a cavity optomechanical sys¬ 
tem with two mechanical modes coupled to the same optical 
cavity mode, where optomechanically-induced transparency 
(OMIT) phenomenon [30-33] occurs. We show that with the 
help of quantum interference, ground state cooling of the me¬ 
chanical mode with uim 'C k can be achieved. Moreover, 
we examine the multiple input cascaded OMIT cooling which 
further suppresses the quantum backaction heating. This ren- 


(a) (b) 



Cooling 


\n+\, 


n + l, 




|n + iSm^,m-l) Heating^ ’'.Control 

Cooling / 

/ - 


l,m + l) 




, [n, fn 


FIG. 1: (color online) (a) Sketch of a typical optomechanical sys¬ 
tem with two mechanical modes b and c coupled to the same optical 
cavity mode a. The cavity is driven by a cooling laser and a control 
laser, (b) Energy level diagram of the system. |n, rric, m) denotes the 
state of n photons, rric c-mode phonons and m 6-mode phonons in 
the displaced frame. The red solid (dashed) arrow denotes the cool¬ 
ing (heating) process of mode 6. The blue dotted arrow denotes the 
control laser enhanced coupling between mechanical mode c and the 
optical cavity mode a. 


ders quantum optomechanics with low optical-Q cavities and 
low mechanical frequency resonators. 


MODEL 

In a generic optomechanical system, as shown in Fig. 1(a), 
we consider an optical cavity mode a coupled to two mechani¬ 
cal resonance modes b and c, where b is the mode to be cooled 
and c is a control mode. The cavity is driven by a cooling laser 
and a control laser, with frequencies ojq and wi, respectively. 
In the frame rotating at the cavity resonance frequency uje, the 
system Hamiltonian reads as 

H = Ht + H,, 

Hb = uJmb^b + ga^a(b + b^} + + H.c.), 

He = tCnicC^c + PcO^aic -F c^) + + H.c.). (1) 
















2 


Here Hb (He) describes the Hamiltonian related with mode b 
(c); Win (Winc) IS the resonance frequency of mode b (c); g and 
gc denote the single-photon optomechanical coupling rates; 
Hq (Oi) represents the driving strength and Aq = Wq — Wc 
(Ai = wi — Wc) is the frequency detuning between the cooling 
(control) laser and the cavity mode. For strong driving, the 
linearized system Hamiltonian is given by 

Hl = LOmb\b^ + [GW4+G(*)*ai](6i + b\) 

+ Winccjci + [G^*^a|-|-G^*^*ai](ci -|- c\). (2) 

Here the operators oi, bi and ci describe the quantum fluc¬ 
tuations around the corresponding classical mean fields af¬ 
ter the linearization; G*^*^ = g{aoe' -f aie and 

Gc*^ = -f are the light-enhanced 

optomechanical coupling strengths, with modified detunings 
Aq = Aq -f Aoni, A'j^ = Ai + Aom and ^om = 2((7^/Wni + 
(7^/winc)( |ao|^ + |ai|^); ao and ai are the intracavity held 
from the contribution of the cooling and control laser inputs; 
K, 7 (= Win/Qm) ^”d 7c (= Wmc/Qnic) *e energy decay 
rates of the modes a, b and c. 


QUANTUM NOISE SPECTRUM 

The optical force acting on mode b takes the form F = 
— [G^*^*ai-|-G*^‘^a|]/a;zpFi where a;zPF is the zero-point fluc¬ 
tuation. The quantum noise spectrum of the optical force 
Sff{^) = f {F{t)F{0)) is calculated to be 

^ — ' 3=0 

^ZPF 

X [« + 5c J2k=o ~ 

where x~^ (w) = H + qI ELo [Xmc(w -f A' - 

K) + X*mci-^ - Xmc (w) = 7cK.th + 

l)|Xmc(w)P -f 7c?^c,th|Xmc (-W) P, Xj^ {^) = 

-z(w -f A')-|-k/ 2 andXmc(w) = -i(w - w^c)-l-7c/2, with 
integers j and k being the summation indices. Here Xj {oj) 
represent the optical response to the input light and Xme (w) 
is the response function of the control mechanical mode; 
nth = — 1] and nc,th = — l] 

are the thermal phonon numbers of modes b and c at the 
environmental temperature T. 

In conventional single mechanical mode approach, the 
quantum noise spectrum exhibits a standard Lorentzian curve 
[15]. However, here due to the interaction between the control 
mechanical and the optical cavity modes, the noise spectrum 
[Eq. (3)] is modified to a non-Lorentzianlineshape. This orig¬ 
inates from the quantum interference manifested by OMIT. As 
shown in Fig. 1(b), the system here contains a series of three- 
level subsystems relevant with OMIT for heating suppression. 
In the presence of the control held, the transition amplitude 
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FIG. 2: (color online) Quantum noise spectra (arbitrary units). The 
solid curves denote (a) S'ff(w), (b) 5 ff(w) and (c) Sff(w) for 
Wm/ F, = 0.02, Wmc/At = 2, Aq= Wm, Aj = Wmc, g /Wm = 10 , 

gcjuJiac = 5 X ao = ai = 10®, Qmc = 10^ and nth = 10®. 
The black dotted curves in (a)-(c) correspond to the results without 
the control mode (gc = 0). In (a) and (b), the position of the sharp 
peaks and dips are marked, (d) Zoom-in view of the shaded region in 
(a)-(c). Spp{uj), blue dashed curve; Spp(cj), green dashed-dotted 
curve; Sfp{oj), red solid curve. The black dotted curves corre¬ 
sponds to the result without the control mode and the control laser 
[S^’f(w)]. The gray vertical lines denote to = ±Wm. 


between the two pathways (red dashed arrow and blue dotted 
arrow) destructively interfere, leading to the suppression of 
the heating transition. 

For the unresolved sideband regime (wm <C k), the quan¬ 
tum noise spectra are plotted in Fig. 2 with parameters 
^m /^ — 0.02, Wmc/Av = 2, Aq = Wm? A^^ = Wmci 
g/Wm = 10“®, 5c/Wmc = 5 X 10“"^, ao = «! = 10®, 
Qmc = 10^ and nth = 10®. Note that Spp(uj) cor¬ 
responds to the spectrum associated with the cooling laser 
(with detuning Aq) and Spp{uj) represents the spectrum 
related to the control laser (with detuning A'^). With¬ 
out the control mechanical mode (gc = 0), they reduce 
to S°p:p{uj) = K|aoXo(w)|^ 5 ^/a;|pF and Sp’p(uj) = 

K \aixi 5^/2^zpF’ which are Forentzians with the cen¬ 
ters at w = —Aq and — A'j^ [black dotted curves in Figs. 2(a) 
and 2(b)], respectively. In contrast, with the presence of the 
control mode, a series of OMIT resonances appear in the spec¬ 
tra. For Spp{uj), those dips/peaks are located at w = ± Wmc, 
—(5± Wmc [Fig. 2(a)], where 5 = Aq — A'^ represents the two- 
photon detuning of the input lasers. Thereinto, the resonances 
at w = ± Wmc originate from the interaction between the cool¬ 
ing laser and mode c, which changes the mode density for 
the cavity held to absorb/emit a phonon with energy Swmd 
the resonances at w = —<5 ± Wmc stem from the interaction 
among the cooling laser, the control laser and mode c. Analo¬ 
gously, for Spp(uj), the resonances are located at w = ± Wmc, 
S ± Wmc [Fig. 2(b)]. For cooling of mode b, the dips/peaks 

























3 



t (units of k'^) 



FIG. 3: (color online) (a) Time evolution of the mean phonon num¬ 
ber nb{t) with control mode (red closed circle) for k = 0.05, 

U.’mc/^ = 2, = U.’m, = CJmc, H/u.’m — 10 , 

= 5 X Qo = 1200, ai = 1600, g^c = lO"*, 
Q-m ~ 10® and nth = 10®. The results without the control mode 
and the control laser (gc = 0, ai = 0) for Aq = —cjm (blue open 
circle) and Aq = tUm (green triangles) are plotted for comparison, 
(b) Same as (a) except that cUm/K = 0.02 and ao = ai = 10®. The 
shaded regions in (a) and (b) denote rib < 1. (c) Cooling rates Topt 
as functions of in the unit of n with the presence (red solid curve) 
and absence (blue dashed curve) of the control mode and the control 
laser; the parameters are the same as Fig. 3(b). 

at w = itwin are relevant, since Azp = Sff (±Wm) a;|pp are 
the rates for absorbing and emitting a &-mode phonon by the 
cavity field, corresponding to the cooling and heating of mode 
b, respectively. With the appropriate value of the two-photon 
detuning at i5 = ujmc + Wm, the OMIT lineshapes in Sff{ijj) 
can be tuned to appear at w = iwm, as shown in Figs. 2(c) 
and 2(d). At a; = —Wm, it exhibits a deep OMIT window, 
which reveals the suppression of heating process, originating 
from the destructive interference. Although a shallow dip also 
appears at a; = Wm, it only slightly decreases the mode den¬ 
sity. The reason is that, with |A'r|» |Aq|, this dip is located 
far away from the center (w = — A'^) of the Lorentzian back¬ 
ground in Spp{u!), as shown in Fig. 2(b). 

COVARIANCE APPROACH 

To verify the destructive quantum interference effect, we 
next solve the quantum master equation and use covariance 
approach [34, 35] to obtain exact numerical results. The mas¬ 
ter equation is given by p = i[p, Hf] + + j{nth + 


l)T>[bi]p+jnth'D[b\]p+jc{nc,th + i)T^[ci]p+'ycnc,th'D[cl]p, 
where Il[o]p = opd^ — {d^dp + pd^6)j2 denotes the the stan¬ 
dard dissipator in Lindblad form. In Fig. 3(a) and 3(b) we 
plot the time evolution of the mean phonon number for 
typical parameters. For single mechanical mode case in the 
unresolved sideband regime, with red detuning input laser, the 
mechanical motion is only slowly cooled with a small cooling 
rate (net optical damping rate) Fopt = A_ — A+, without 
reaching the ground state. However, in the presence of the 
control mode and the control laser, the cooling rate can be en¬ 
hanced for more than two orders of magnitude [Fig. 3(c)], 
and ground state cooling with mean phonon number nf, < 1 
is achievable, even for sideband resolution Wni/w as small as 
0.02. It should be emphasized that in this case the cooling 
laser is blue detuned, which is quite different from the single 
mechanical mode approach. For the latter, blue detuning leads 
to amplification instead of cooling of the mechanical motion. 
The blue-detuning cooling is the unique property originating 
from the quantum interference which modifies the noise spec¬ 
trum of the optical force. The OMIT lineshape can be viewed 
as the inverse of the standard Lorentzian, thus the detunings 
for cooling are just opposite to that for the single mechanical 
mode case. 

CASCADED OMIT COOLING 

To further suppress quantum backaction heating, we pro¬ 
pose the use of additional coherent laser inputs, resulting in 
cascaded OMIT cooling. For N inputs, the quantum noise 
spectrum of the optical force takes the same form as Eq. (3) 
except that the summation indices (j, k) run from 0 to A — 1. 
As displayed in Figs. 4(a)-4(d), for two inputs, the suppres¬ 
sion of heating for mode b is the contribution of suppressed 
= Spp ( x^pp, while = Spp (u^m)2^zPF 
also slightly suppressed. In the presence of the third input 
laser with detuning A 2 = Aq — 2(wnic AWm), the interac¬ 
tion involved with the control mode results in the suppression 
of A_|_ = Spp{ and A_ = Sppi^uj^^x^pp. This 

results in the net optical damping rate Fopt = Ylk=oi^- ~ 
Ap) ~ A°_ — Ap [Fig. 4(c) and (d)]. More generally, for 
N inputs with detuning A], = Ag — kitOrac + i^m) (fc = 
1,2, ...A— 1), we obtain Fopt — A^_—Ap~"^. Note that the re¬ 
maining heating rate A^~^ is much smaller than the original 
heating rate due to the large detuning AJ.for the (A — 1 )- 
th input. In Fig. 4(e) we compare the cooling dynamics 
between two inputs and three inputs for typical parameters, 
which shows that the cascaded OMIT cooling enables larger 
cooling rate and lower cooling limit, with ground state cooling 
achievable even for k, = 0.01. 

COOLING LIMITS 

In Fig. 5 the fundamental cooling limits riniin as func¬ 
tions of the sideband resolution k are plotted. The exact 
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FIG. 4: (color online) (a)-(b): Quantum noise spectra (arbitrary 
units) for two inputs (a) and three inputs (b). 5 'ff(w), red solid 
curve; Spp{uj), blue dashed curve; Spp{uj), green dotted curve; 
Spp{u!), purple dashed-dotted curve. The gray vertical lines denote 
u = itUm. (c) Scattering interpretation of optomechanical interac¬ 
tions with three inputs. The shaded Lorentzian represents the mode 
density of the optical cavity. The green vertical arrows denote the in¬ 
put lasers. The solid (dashed) curved arrows denote the anti-Stokes 
(Stokes) scattering processes relevant with mode b (red) and mode c 
(blue), (d) Energy levels and interpretation of the suppression (de¬ 
noted by the “x”) of heating and cooling with three inputs, (e) Time 
evolution of the mean phonon number nb{t) with two inputs (blue 
open circle) and three inputs (red closed circle). The shaded region 
denotes rib < 1- Here 0.-2 = 10®, cumc/K = 2, = 0.01 and 

other parameters are the same as Fig. 3 (b). 


FIG. 5; (color online) Fundamental cooling limits rtmin as functions 
of uJhi/k for two inputs (red closed circles) and three inputs (blue 
open circles). The result for single mechanical mode approach (black 
dotted curve) is plotted for comparison. The shaded region denotes 
n-rnin < 1. Here ao = Qi = 0:2 = 500 and other parameters are the 
same as Fig. 4. 

EXPERIMENTAL EEASIBILITY 

It should be stressed that the OMIT cooling described here 
adds little complexity to the existing optomechanical system, 
which is crucial in the experimental point of view. Compared 
with the conventional backaction cooling approach, the ad¬ 
ditional requirement here is a control mechanical mode and 
one (or more) input laser. It is experimentally feasible for 
various optomechanical systems within current technical con¬ 
ditions. On one hand, many optomechanical systems pos¬ 
sess abundant mechanical modes with different resonance fre¬ 
quencies, since the oscillation have different types and orders. 
This situation can be found in optomechanical systems us¬ 
ing whispering-gallery microcavities [36, 37], photonic crys¬ 
tal cavities [38, 39], membranes [40,41], nanostrings [42] and 
nanorods [43, 44] amongst others. Usually only one mechan¬ 
ical mode is used in most optomechanical experiments, while 
exciting an additional mechanical mode is often unintended. 
On the other hand, composite optomechanical systems, con¬ 
taining two independent mechanical resonators, are also con¬ 
ceivable. For example, in Fabry-Perot cavities, the motion 
of one mirror acts as an control mechanical mode while the 
other mirror is to be cooled [Fig. 1(a)]. In the near-held op¬ 
tomechanical system [42], to cool the nanostrings, the control 
mode can be selected from the vibration of the microtoroid. 


numerical results are obtained from master equation simula¬ 
tions. The black dotted curve shows the best result for con¬ 
ventional single mechanical mode approach, given by riniin = 
At/(Twm), which is obtained when Aq = — k/ 2[14]. Itreveals 
the great advantage of OMIT cooling and cascaded OMIT 
cooling, with possibility for ground state cooling even when 
Wm/tt ~ 3 X 10”^, which goes beyond the resolved sideband 
limit by nearly 3 orders of magnitude. Note that Fig. 5 shows 
the cooling limits increase as Wm/tt increases from ~ 0.02 
to a larger value. This is a result of blue detuning induced 
heating, which becomes significant when Wm/At is large. 


CONCLUSION 

In summary, we have presented the OMIT cooling scheme 
allowing ground state cooling of mechanical resonators be¬ 
yond the resolved sideband limit. It is demonstrated that by 
employing the OMIT interference, quantum backaction heat¬ 
ing can be largely suppressed, extending the fundamental limit 
of backaction cooling. The scheme is experimentally feasible, 
which requires another control mechanical mode and multiple 
laser inputs. Such a M-O-M system (M, mechanical mode; O, 
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optical mode) studied here offers potential for cooling enor¬ 
mous mass scale resonators [45, 46], which possess small 
resonance frequencies. Together with the recently examined 
multi-optical-mode [28, 47-50] and multi-mechanical-mode 
[51-55] systems,it is shown that such interference effect in 
multi-mode cavity optomechanics provides unique advantage 
for both fundamental studies and broad applications. Recently 
we noticed a related work [56], but here we use the covari¬ 
ance approach to examine the fundamental cooling limits and 
present detailed analysis of cascaded OMIT cooling. This 
paves the way for the manipulation of macroscopic mechani¬ 
cal resonators in the quantum regime. 
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